Long-range attraction between probe particles mediated by a driven fluid 
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The effective interaction between two probe particles in a one-dimensional driven system is studied. 
The analysis is carried out using an asymmetric simple exclusion process with nearest-neighbor 
interactions. It is found that the driven fluid mediates an effective long-range attraction between 
the two probes, with a force that decays at large distances x as —b/x, where b is a function of the 
interaction parameters. Depending on the amplitude b the two probes may form one of three states: 
(a) an unbound state, where the distance grows diffusively with time; (b) a weakly bound state, in 
which the distance grows sub-diffusively; and (c) a strongly bound state, where the average distance 
stays finite in the long time limit. Similar results are found for the behavior of any finite number of 
probes. 

PACS numbers: 02.50.Ey, 05.40.-a 



Probe particles are a powerful tool in the study of 
properties of solutions. Interactions between solute par- 
ticle arise from direct interactions as well as a solvent- 
mediated part to which the free energy of solvation con- 
tributes. The latter contribution is sometimes rather 
significant, with surprising effects such as the hydropho- 
bic attraction [l| in which entropy plays the dominant 
role. Strong entropic effects have also been observed in 
anisotropic fluids |2j. In equilibrium one frequently ap- 
proximates these interactions by setting up a suitable 
model and calculating an effective potential of mean force 
between solute particles, using the potential distribu- 
tion theorem [2|. While more quantitative studies re- 
quire quite sophisticated modelling, simple lattice mod- 
els often suffice to successfully explain generic features of 
these interactions. Difficulties of conceptual rather than 
practical nature, however, arise in equilibrium-based ap- 
proximation schemes if the solvent is in a strongly non- 
equilibrium state and therefore no notion of free energy 
exists. Analyzing the interactions between probe par- 
ticles in non-equilibrium fluids would thus be of great 
interest. 

Far from equilibrium driven systems of particles mov- 
ing steadily have been a subject of extensive studies in 
recent years Q . The minimal model which has been used 
to describe these systems is the asymmetric simple ex- 
clusion process, whereby self-avoiding particles hop on a 
lattice with rates which favor motion along the direction 
of the drive [MEl- In one dimension (Id) these studies re- 
sulted in detailed calculations of a variety of steady state 
properties, including phase diagrams, density correlation 
functions and other collective features. In particular, it 
has been observed that a probe particle introduced into a 
Id driven fluid is attracted to regions of large density gra- 
dient (on molecular scale) in the steady state densitypro- 
file, thus serving as a microscopic marker of shocks |jj . A 
local shock is accompanied by a long-range algebraically- 
decaying density profile away from the probe. When the 



system contains two probe particles, the density profile 
induces a long-range attractive interaction between the 
probes. For example, in the case of the totally asymmet- 
ric simple exclusion process (TASEP) on a ring with two 
'second-class' probe particles, the steady-state distribu- 
tion of the distance x between the two probes was found 
to decay as x~ 3 / 2 for large x @|. This implies that the 
two probes form a weakly bound state, where the average 
distance is infinite. This exactly soluble case is very spe- 
cific though, as the second-class particles do not influence 
the motion of the fluid particles. On the other hand, it 
is known that the direct interaction between solute and 
solvent could have a strong dynamical effects, and should 
not in general be ignored. For example, direct interaction 
between solute and solvent may lead to a non-monotonic 
dependence of solute and solvent diffusion coefficients on 
the solute diameter 

In this Letter we present a dynamical approach to the 
problem of probe particles in a Id driven fluid. This ap- 
proach is applicable to a broad class of Id systems, which 
in addition to excluded volume display short range inter- 
actions, and where probes directly influence the motion 
of the fluid particles. Moreover, it enables us to study 
not only the steady state properties, but also the tempo- 
ral relaxation to the steady state. We find that within 
this broader class of systems, two probe particles may 
form one of three states: unbound, weakly bound, or a 
strongly bound state. In the unbound state the average 
distance between the probe particles grows diffusively as 
y/i, as is the case for non-interacting particles. In the 
weakly bound state, the steady-state distance distribu- 
tion decays algebraically for large x with a power-law 
x~ a , where the exponent 1 < a < 2 is a function of the in- 
teraction parameters. The steady state average distance 
is infinite. The approach to steady state is sub-diffusive, 
whereby the average distance grows as t u with v < 1/2. 
In the strongly bound state the average distance decays 
algebraically in time to a finite value at steady state. In 
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the steady-state, the distance distribution takes the form 
x~ a , with a > 2. The case of more than two probe par- 
ticles may also be studied within the approach suggested 
in this Letter. 

We now introduce a model within which the interaction 
between probes may be analyzed. The model is defined 
on a Id lattice of ring geometry, where each site can 
be occupied by either a positive particle (+), a negative 
particle (— ), or it may be occupied by a probe particle 
(0). Apart from the exclusion interaction, the 'charged' 
particles are subject to a nearest-neighbor interaction, 
defined by the potential 
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where s$ = 0,±1 according to the occupation of site i, 
and — 0.8<e<lis the coupling constant fl(ij . The 
dynamics of the model is defined by a random-sequential 
local dynamics, whereby a pair of nearest-neighbor sites 
is selected at random, and the particles are exchanged 
with rates 
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with rate 1 + AV 

with rate 1 . (2) 



Here AV is the difference in the potential V between 
the initial and final states. This dynamics conserves the 
number of particles of each species, and is symmetric 
under the exchange of charges and spatial direction. The 
model is a generalization of , which has been studied 
in 0, 0] . In the case of no nearest neighbor interaction, 
e = 0, the dynamics defined above reduces to that of the 
TASEP with second-class particles 0, II] which are the 
probes. 

In the case of no probe particles, the steady state of 
this model has an Ising-like measure 01 • This is also 
expected to be the local steady-state measure away from 
any probe when the density of probes is zero. We consider 
first the case of two probe particles in the system. Our 
aim is to derive dynamical equations for the two probes, 
and study the evolution in time of the distance n between 
them, starting from n = 0. 

For e = it has been shown [9J that the statistical 
weight of all configurations with a given n is proportional 
to Z n , the partition function of the TASEP on an open 
chain of length n in the maximal current phase 5] . This 
observation has been used to show that the probability 
to find the two probes at a distance n from each other 
decays for large n as n~ 3 / 2 8]. In addition, this result 
can be used to estimate the currents of particles which go 
in and out of the segment trapped between the two probe 
particles. One finds that the outgoing current of + (— ) 
particles through the right (left) probe takes the form of 
the steady state current of the TASEP, given for large n 
by jn Ut = i fl + Hr) • ^hc PP 0Sm g currents, namely 
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FIG. 1: (a) Schematic diagram of the positive (black) and 
negative (gray) currents, as seen from the two probe particles. 
Here n is the distance between the two probes. For e = one 
has j = 1/4 and b = 3/2. (b) A random-walk representa- 
tion of the dynamics of the two probe particles, as deduced 
from the currents in (a). Each probe carries out forward and 
backward hops with the indicated rates. 



that of — (+) particles incoming through the right (left) 



probe, take the form j ln = \ ( 1 
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N _ 2 I ■ in the limit 
we are concerned with, namely N 3> n and N — > oo, this 
current is well approximated by j ln — j [Fig.^a)]. 

These considerations have been extended to the case 
e ^ 12]. Following [13] one expects that the current of 
+ (— ) particles bypassing the right (left) probe takes the 
same form as that of the current in an open segment of 
the same length, governed by the same dynamics. This 
current is given by j° ut = jo(e) (1 + b(e)/n), where 
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1. For the relevant values of e one has 

< b < 9/4. Also, similar to the e = case, the incoming 
current is given in the large N limit by f n = Jo(e). 

The considerations described above may be used to 
derive dynamical equations for the two probes. These 
considerations suggest that the two probe particles be- 
have as two coupled random walkers. Each probe moves 
away from the other probe with rate j Q , and move to- 
wards it with rate j (1 + b/n) [Fig.^b)]. Thus, one can 
model the time evolution of the distance n between the 
probes by the master equation 



dP(n, t) 



dt 



P(n - l,t) - P(n,t) 



1 + — b —)p(n+l,t)- (1 + - ) P(n,t) (4a) 
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for n > 0, and 



dP(0,t) 



= -P(0,t) + (l + b)P(l,t) 



(4b) 



for n — 0. Here the time t is a rescaled by a factor 2j . 
For large n one can use the continuum limit, and can 
be rewritten as a Fokker-Planck equation (FPE) 
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with the boundary conditions P(x,t) — > as a; — ► cxd. 
It is easy to see that the steady-state solution of this 
equation is of the form 



(6) 



This in agreement with the exact result x~ 3 / 2 in the case 
e = o||. Note that the force term b/x in the FPE JSJ) cor- 
responds to a logarithmically increasing fluid-mediated 
effective potential U ~ In a; between the probe particles. 

In order to study the approach to steady state we 
make the scaling hypothesis P(x,t) — x t~" f(x/t a ), 
with f(u) a scaling function, which satisfies the bound- 
ary condition /(0) = const., and the normalization con- 
dition Jg 00 P(x) Ax — 1. Here a cutoff 9 is introduced to 
prevent divergence at x — ► for b > 1. Substituting this 
form into JSJ) one finds that a non-trivial solution exists 
only when a = 1/2. In the range b > 1 the normalization 
condition yields (3 — 0. Defining u = x/^/t, the scaling 
function f(u) satisfies for b > 1 
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b)f(u) + uf"(u) = 0. 
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This equation reduces to f'(u) ~ u h exp(— it 2 /4), and its 
solution is given by f(u) ~ T^^T'l")' wnere Tl !^) 
is the incomplete T-function. In the range b < 1 the 
normalization condition yields /3 = (1 — b)/2 and hence 



) + «/» = - 



1-6 



''/(«) (8) 



which is solved by /(«) ~ exp(— u 2 /4). 

We now turn to calculate the average distance between 
the particles, (x(t)) = dx xP(x,t). Inspecting the 
large t behavior of this integral one finds 



s(*)> 



f t 1 ' 2 b < 1 
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t 1 - 6 / 2 1< b < 2 

log(t) 6 = 2 

J 4 + Si-( fc / 2 - 1 ) 6>2 
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Here A and -B are some non-universal constants. It is ev- 
ident that depending on the value of b three regimes can 
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FIG. 2: Results of numerical simulations for the average dis- 
tance between the two probes for several values of e, and thus 
of b. Lines are drawn with the expected slope, according to 
Eq.0 



be identified. For b < 1 the two probes behave as two de- 
coupled particles, and the distance between them grows 
diffusively. For 1 < b < 2 the probes are weakly bound, 
and the distance between them diverges sub-diffusively. 
For b > 2 the two probes are strongly bound, and the 
average distance between them approaches a constant al- 
gebraically in time. 

To test the validity of these results we carried out nu- 
merical simulations of the model with two probe parti- 
cles. The simulations were performed on a ring of length 
N = 10 5 , with equal number of + and — particles. Start- 
ing with the two probes located on nearest neighbor sites, 
we followed the evolution of their distance n(t). The av- 
erage distance (n(t)) was obtained by averaging over 5000 
realizations of the noise. Results for different values of b 
in the three regimes are presented in Fig. [2 These results 
agree very well with the theoretical predictions. 

We now consider the case of a finite number M of 
probes, and show that the scaling behavior, Eq.|5| holds 
for the average distance between any pair of probes. The 
evolution of M probes can be be modelled as a zero- 
range process (ZRP) 01, defined on a one-dimensional 
ring of M sites, where each site can be occupied by any 
number of particles. At any time step a site is chosen 
at random, and one of its particles hops to either one 
of its two nearest neighbors with the same rate w n — 
jo (1 + b/n). Here n is the occupation of the departure 
site. To make the correspondence between this ZRP and 
the driven model (0 we identify the occupation number 
of site i with the distance between the probes i and i + 1 
in the driven model. Moreover, since we are interested 
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in the limit N — > oo, we take site i = M to be occupied 
by an infinite number of particles. This site serves as a 
reservoir, omitting particles with current j to both sides. 
For this zero-range process it can be shown 0] that the 
distribution of occupation numbers in the steady state 
takes the form of a product measure, 

Af-l 

P s t {n-i,n 2 , ■ ■ ■ ,n M -i) ~ JJ nT b . (10) 

i=l 

In the particle system, this steady-state distribution 
holds exactly only in the case e = 0. However, it has 
been shown that the correlations between adjacent seg- 
ments can be ignored at large distances m also for e ^ 
0)0]. With the steady-state distribution at hand, the 
temporal approach to steady-state can be studied in a 
similar manner to the case of two probes. In particular, 
we make the scaling hypothesis 

P(xi,-- ■ ,x M -i,t) = 

^-^-(^/(l-,^) , (ID 

where, again, Xi is the continuous variable correspond- 
ing to rii. A straightforward analysis of this scaling form 
shows that holds for all (x^ , and thus for the distance 
between any pair of probes. This agrees well with numer- 
ical simulations of systems with 5 probes, supporting the 
scaling hypothesis (details are not presented). 

Our discussion so far was limited to the case where 
the number of positive particles was equal the number 
of negative particles. In it has been argued that 
the physical picture of domain currents underlying the 
dynamical approach developed above still holds when the 
densities of the two species are not equal. In this case, the 
amplitude b becomes a function of both e and rj, where 
r] is the density of, say, the positive particles (see for 
explicit expressions). It is evident that b{e,rj) increases 
when rj deviates from |. Thus, increasing the density of 
one species at the expense of the other results in larger 
6, and may change the state of the probes from unbound 
to weakly bound to strongly bound state. 

In summary, we introduced in this Letter a dynam- 
ical approach for the study of probe particles in one- 
dimensional driven fluids. This approach generalizes the 
exact result of §| , obtained for the steady-state distribu- 
tion of two probe particles which do not affect the motion 
of the fluid, to a broader class of models, which are not 
exactly solvable. Here any finite number of probe par- 
ticles, which influence the motion of the fluid, are con- 
sidered. Moreover, the dynamical approach enables one 
to study not only steady-state properties, but also the 
long-time temporal behavior of the probes. It is found 
that a logarithmically increasing fluid-mediated effective 
potential U ~ In a; acts between the probes when they 
are a distance x apart. It is remarkable that the probe 



particles behave as if they were in thermal equilibrium 
even though the fluid is very far from equilibrium. This 
is reminiscent of the motion of a boundary-induced shock 
in non-conserving driven systems in contact with a reser- 
voir. The position of the shock has recently been found 
to perform Brownian motion in an effective potential dy- 
namically generated by the fluid medium 16]. 
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